Time dependence of joint entropy of oscillating quantum systems 
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Abstract 

The time dependent entropy (or Leipnik's entropy) of harmonic and damped harmonic oscillators 
is extensively investigated by using time dependent wave function obtained by the Feynman path 
integral method. Our results for simple harmonic oscillator are in agrement with the literature. 
However, the joint entropy of damped harmonic oscillator shows remarkable discontinuity with time 
for certain values of damping factor. According to the results, the envelop of the joint entropy 
curve increases with time monotonically. This results is the general properties of the envelop of 
the joint entropy curve for quantum systems. 
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I. INTRODUCTION 



The investigation of time dependent entropy of the quantum mechanical systems attracts 
much attention in recent years. For both open and closed quantum systems, the different 
information-theoretic entropy measures have been discussed [2j, y, m- In contrast, the 

nn 

joint entropy [5|, [6|] can also be used to measure the loss of information, related to evolving 
pure quantum states The joint entropy of the physical systems which are named MACS 
(maximal classical states) were conjectured by Dunkel and Trigger 8]. According to Ref. Q], 
the joint entropy of the quantum mechanical systems increase monotonically with time but 
this results are not sufficient for simple harmonic oscillator 

The aim of this study is to calculate the complete joint entropy information analytically 
for simple harmonic and damped harmonic oscillator systems. 

This paper is organized as follows. In section II, we explain fundamental definitions 
needed for the calculations. In section III, we deal with calculation and results for harmonic 
oscillator systems. Moreover, we obtain the analytical solution of Kernel, wave function in 
both coordinate and momentum space and its joint entropy. We also obtain same quantities 
for damped harmonic oscillator case. Finally, we present the conclusion in section IV. 

II. FUNDAMENTAL DEFINITIONS 

We deal with a classical system with d = sN degrees of freedom, where N is the particle 
number and s is number of spatial dimensions 8|. We assume that the density function 
g(x,p,t) = g(xi, Xd,Pi, ■■■,Pd, t) which is the non-negative time dependent phase space 
density function of the system has been normalized to unity, 



dxdpg(x,p,t) = 1. (1) 

The Gibbs-Shannon entropy is described by 

S(t) = -^y J dxdpg(x,p,t)ln(h d g(x,p,t)), (2) 

where h = 2irh is the Planck constant. Schrodinger wave equation with the Born interpre- 
tation 10j is given by 



in^ = m. (3) 
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The quantum probability densities are denned in position and momentum spaces as 
and \ip(p, t)\ 2 , where \ip(p, t)\ 2 is given as 

dxe~ ipx/h 



Leipnik proposed the product function as 



(27rh) d /' 



-ip(x,t). 



(4) 



gj (x, P ,t) = \ij(x,t)\ 2 \^p,t)\ 2 >0. 



(5) 



Substituting Eq. (5) into Eq. (2), we get the joint entropy Sj(t) for the pure state ip(x,t) 
or equivalently it can be written in the following form [8] 



(6) 



Sfi) = - J dx\i,(x,t)\ 2 \n\i,(x,t)\*- I dp\^(p,t)\ 2 ln\^(p,t)\^ 
- \iyh d . 



We find time dependent wave function by means of the Feynman path integral which has 

n 

form [11] 

px"=x(t") 

K(x",t";x',t') = / DxtyeWH 

Jx'=x(t>) 



Dx(t)eKb L ^ dt . (7) 
The Feynman kernel can be related to the time dependent Schrodinger's wave function 

oo 

K(x", t"; x', t>) = J2 r n (x', t')Mx", t"). (8) 



n=0 



The propagator in semiclassical approximation reads 

i d 2 



K(x", t"; x', t') 



2nh dx'dx" 



l 1/2 e iS cl (x",t";x',f)_ 



(9) 



The prefactor is often referred to as the Van Vleck-Pauli-Morette determinant 
F(x", t"; x', t') is given by 
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13). The 



F(x", t"; x', t') 



i d 2 

2nh dx'dx" 



Sci{% -,t j x , t ^ 



1/2 



(10) 
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III. CALCULATION AND RESULTS 



A. Simple Harmonic Oscillator (SHO) 

To get the path integral solution for the SHO, we must calculate its action function. The 
Lagrangian of the system is given by 

L(x,x,t) = ^(x 2 - 1 -u 2 x 2 ) (11) 

Following a straightforward calculation, it is given by: 

S(x d (n x cl (t')) = ^-\{x'2 + x' 2 ) coscot - 2x' d x%\ (12) 
z sin uiz 

with t = t" — V and x' d = x ,x d = x. Substituting Eq. (9) into Eq. (7), we obtain the 



Feynman kernel 



muj . i , mw., 9 9 . 2xnX . „. 

K(x, x 0] t) = — — 2 exp{-— x 2 + x 2 ) cot out - }. 13 

27rm sin ut 2tn sm cjt 

By the use of the Mehler-formula 

00 " " ' „2 , „,2Wi , ^2 



e -(x 2 +?, 2 )/2 



El ,2,n , \ tt i \ 1 r 4a;wz - (x 2 + y)(l + z ), 
— (-) 2 HJx)H n (y) = exp — ^ ^ (14) 

nP2 ; v 7 vy; vT^2 P[ 2(1 -z 2 ) J v ; 



n=0 



where if n is Hermite polynomials, we can write the Feynman kernel defining x = ^mu/hxo, 
y = a/ mu j fix and z = e~ lwT 

oo 

AT(x, x ; t) = J2 e- itEn/h y*{x )y(x) (15) 

n=0 

with energy-spectrum and wave-functions: 

E n = hu(n+^), (16) 



, . . muj , i . muo . , muj 9x . . 

Time dependent wave function of the SHO is defined as 

m(x,t) = J K(x,x ;t)^(x ,0)dx . (18) 

It can be written as 

/muj\ 1 / 4 r a a 2 iut} r a 2 9 ,- . ... 

= (^rj ex P { - 4 - y - - g-} ex P [ - x e + ac * e J (19) 



where x or a is mean of the Gaussian curve. The probability density has 

"muj \ V 2 



l*(M)l a = f 



cxp 



(a — a cos uty 



(20) 



where a = \P^rX. Thus it can be written as 

V n 



|*(x,t)| J 



71"^ 



cxp 



X — X COSLdty 



(21) 



This has been shown in FigJTJ In momentum space, the probability density has the form 



\i>(p,t)f 



(-L- 

\mujT:Ti 



1/2 



exp 



+ 



mux 2 



.rntuh 2h 
The joint entropy of harmonic oscillator becomes 

Amu 



cos 2a; (t) 



2px 



smuj(t) 



(22) 



Sj(t) 



In - + 

2 h 



x 2 sin 2 uj{t). 



(23) 



In FigJ2l the joint entropy of this system was plotted by using Mathematica in three dimen- 
sion. As known from fundamental quantum mechanics and classical dynamics, displacement 
of simple harmonic oscillator from equilibrium depends on harmonic functions (e.g sine or 
cosine function). Therefore, other properties of the SHO systems indicate the same har- 
monic behavior. If the frequency of the SHO is sufficiently small, the system shows the 
same behavior as the free particle 8| . As seen from Figj3] and FigJH envelop of the sinusoidal 
curve is also monotonically increase with omega and constant with time at constant omega, 
respectively. When the frequency increases, the joint entropy of this system indicates a fluc- 
tuation with increasing amplitude with time. If t goes to zero, it is important that Eq.(20) 
is in agreement with following general inequality for the joint entropy: 

Sj(t) > ln(|) (24) 

originally derived by Leipnik for arbitrary one-dimensional one-particle wave functions. 



B. Damped Harmonic Oscillator (DHO) 



The DHO is very important physical system in all physical systems defining an interaction 
with its environment. The Lagrangian of the DHO is given by 



L{x, x, t) = (^x 2 - ^u 2 x 2 + j(t)xj) . 



(25) 
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The time dependent wave function is obtained as 



ip n {x,t) = I cLxqK{x 1 t] xq, 0)i/j(x, 0) 



N 



(2 n n!) l/2 eX Pi ' 



n + - ) cot 1 x 
2/ 



x ^- + cotut + fj jexp[-(Ac 2 + 
+ 2Bx)]H n [D(x - E)}. 



(38) 



To simplify the evaluation, we set j(t) = 0. Such that kernel and wave function of the DHO 
become 



K(x,t;x ,0) = 



( mue^ 2 \V2 



exp 



vm 



7 (x 2 - A 2 ) + - 
4n V sinwt 



2oo 



2irih sin cut 
x [(x^ + x 2 e 7 *)coscjt-2e 7i/2 xx ] 



x 



(39) 



where u = (u 2 — 7 2 /4) 1 / 2 and 

, , s N 
1p n {X,t) 



( 2 n n !)l/2 

Where D, A and N are 



exp < — % 



n H — ^ cot 1 f — + cot tut 
2/ \2w 



D(t) 



ae 



it 12 



r](t) sin cut ' 



v 2 (t) 



7 7 2 
H — cos cut + esc cot, 



Au 2 u 



A(t) 



muo 



2h 



?7 2 (t)sin uot \2uj 



7 ■y/2u + cotuot 
+ * ( cot uot H , 2 



tf n [L>a;]exp[-Ar 2 ]. (40) 

(41) 
(42) 

(43) 



r? 2 sin cut 



and 



muo\ 1 / 4 exp (2^ 



7T^, / ^(t) (sin cut) 1 / 2 
The ground state wave function is given by 

ipo(x,t) = iV(t)exp{ -i cot -1 + cot cut) j exp[-A(t) 

So the probability distribution in coordinate space becomes 



\Mx,t)\ = N(ty exp[-2 A' (t)x 2 } 



where A' is defined by 



A\t) 



muo 



2k l7] 2 (t)sm 2 ut 



(44) 



(45) 



(46) 



(47) 
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The probability density in coordinate space is shown in Figj5] and Figj6] for the different 
values of 7. The probability density in momentum space can be written easily 

The time dependent joint entropy can be obtained from Eq. (2) as 



m - w?^[q*w? - \) - IJumw^mm ~ s)J - " lM49) 

The joint entropy depends on damping factor 7. When 7^0, all the above results 
are converged to simple harmonic oscillator. However, when the 7 7^ 0, the joint entropy 
has remarkably different features of the SHO. As can be seen in FigJ7]and FigJHl the joint 
entropy of the DHO has very interesting properties. One of the most important properties 
of the joint entropy is the probability of taking values for small 7 values. As we know 
from literature the joint entropy must be positive and monotonically increase. However, this 
system has different properties from literature because of periodically discontinuity of the 
joint entropy. On the other hand, envelop of this curve is also monotonically increase with 
time for large 7. As can be shown these results, the envelop of the joint entropy curves 
has general properties as monotonically increase for quantum systems. Thus, we have found 
that the joint entropy is depend on properties of investigated system. 



IV. CONCLUSION 



We have investigated the joint entropy for explicit time dependent solution of one- 
dimensional harmonic oscillators. We have obtained the time dependent wave function by 
means of Feynmann Path integral technique. Our results show that in the simple harmonic 
oscillator case, the joint entropy fluctuated with time and frequency. This result indicates 
that the information periodically transfer between harmonic oscillators. 

On the other hand, in the DHO case, the joint entropy shows a remarkable smooth 
discontinuities with time. It also depends on choice of initial values of parameter i.e. u. 
These results can be explained as the information exchange between harmonic oscillator 
and system which is supplied damping. But the information exchange appears in certain 
values of time for damping. If the damping factor increases, the information entropy has 
not periodicity anymore. Moreover, for certain values of the damping factor, the transfer of 
information between systems is exhausted. 
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FIG. 1: \^(x,t)\ 2 versus time and coordinate 
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FIG. 3: The joint entropy of simple harmonic oscillator versus uo . 
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FIG. 4: The joint entropy of simple harmonic oscillator versus time. 
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FIG. 5: The probability function as a function of time and coordinate at 7 = 0.1. 
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FIG. 6: The probability function as a function of time and coordinate at 7 = 0.5. 
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FIG. 7: The 3D graph of the joint entropy of damped harmonic oscillator for damping factor^) 
at ujq = 2. 
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FIG. 8: The 3D graph of the joint entropy of damped harmonic oscillator for damping factor^) 
at ujq = 1. 
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